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Abstract. We study the neutrino emissivity of strongly 
magnetized neutron stars due to the charged and neu- 
tral current couplings of neutrinos to baryons in strong 
magnetic fields. The leading order neutral current pro- 
cess is the one-body neutrino-pair bremsstrahlung, which 
does not have an analogue in the zero field limit. The 
leading order charged current reaction is the known gen- 
eralization of the direct Urea processes to strong mag- 
netic fields. While for superstrong magnetic fields in ex- 
cess of 10 18 G the direct Urea process dominates the one- 
body bremsstrahlung, we find that for fields on the order 
10 16 — 10 17 G and temperatures a few times 10 9 K the one- 
body bremsstrahlung is the dominant process. Numerical 
computation of the resulting emissivity, based on a simple 
parametrization of the equation of state of the npe-matter 
in a strong magnetic field, shows that the emissivity of this 
reaction is of the same order of magnitude as that of the 
modified Urea process in the zero field limit. 
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1. Introduction 

It is now well established that the neutron stars which 
are observed as radio-pulsars posses B-fields of the order 
of 10 12 -10 13 G at the surface. The interior fields are un- 
known, but can be by several orders of magnitude larger 
than the ones inferred for the surface. The scalar virial the- 
orem sets an upper limit on the magnetic field stre ngth o f 
a neutron star of the order of 10 18 G (Lai & Shapiro 1991 ). 
Similar conclusion is reached through more sophisticated 
numerical studies (Bocquet et al. 1995 ). 

Recent measurements of the spin-down timescales of 
several soft gamma-ray repeaters, such as SGR 0526-66 
(Mazets et al. |1979| ), SGR 1806-20 (Mu rakam i et al. |1994[ ), 
and SGR 1900+14 (Kouveliotou et al. [1998D with RXTE, 
ASCA and BeppoSAX have made a strong case for SGRs 
as being newly born neutron stars that have very large 
surface magnetic fields (up to 10 15 G). The subsequent 
discoveries of the SGR 1627-41 by BATSE (Woods et al. 



1999|) and SGR 18 01-23 by Ulysses, BATSE, and KONUS- 
Wind (Cline et al. |l999f ) lent further support to the iden- 
tification of SGRs with highly magnetized neutron stars. 
These objects were naturally related to the magnetars, 
which are thought to be remnants of a supernova explosion 
which develop high magnetic fields via a dynamo mech- 
anis m (D uncan & Thompson 1992] , Thompson & Dun- 



1995 ). The magnetars also serve as a model for the 
anomalous X-ray pulsars (AXP) (van Paradijs, Taam & 



van den Heuvel 
thelf, Vasisht, I 
et al. 



1995) such as IE 1841-045 (Kes 73) (Got- 



Dotani 



1999), RX J0720.4-3125 (Haberl 



1997), and IE 2259+586 (Rho & Petre 1997) 



Neutrino-nucleon interactions in the strong magnetic 
fields have been studied recently both in the supernova 
and neutron star contexts. It has been pointed out that 
the neutrino emission from proto-neutron stars, which is 
anisotropic in strong B-fields, could produce the "pulsar 
kicks" if the fields are in excess of 10 16 G (Horowitz & 
Li 199? ; Arras & Lai 199£ and references therein). The 
strong magnetic fields relax the kinematical constrains on 
the direct Urea process and hence give rise to finite neu- 
trino emissivity even when the proton fraction is small 
(Leinson & Perez 1998, Bandyopadhyay et al. 1998] ). The 
effect is most pronounced in the ultra-high magnetic fields 
when protons and the electrons occupy the lowest Lan- 
dau levels. The direct Urea process for arbitrary magnetic 
fields, when the protons and electrons are allowed to oc- 
cupy many Landau levels, has also been studied (Baiko & 
Yakovlev [l999| ). 



The neutrino emissivities via the one-body processes 
sensitively depend on the abundances of the various 
species of baryons and leptons which are controlled by 
the equation of state (EoS) of the dense matter in strong 
magnetic fields. The strong magnetic fields lead to an in- 
crease of the proton fraction (Broderick, Prakash & Lat- 
timer 2000). The muon production and pion condensation 
in strong magnetic softens the EoS of the dense matter 
(Suh & Mathews |l999| ). For the purpose of estimating the 
magnitude of the neutrino emissivities, we employ in this 
paper a simple parametrization of the EoS for the npe- 
matter in strong magnetic fields. 
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The main objective of this paper is to show that the 
strong magnetic fields open a new channel of neutrino 
emission via one-body neutral current brcmsstrahlung, 
which does not have an analogue in the zero field limit. 
We also briefly discuss the direct Urea process, which is 
forbidden in the low-density zero-field limit (as long as 
the triangular condition PFp + PFe > PFn is not sat- 
isfied), but is allowed in strong magnetic fields because 
of the relaxation of the kinematical constrains. We com- 
pare the emissivities of various reactions using a simple 
parametrization of the EoS for the npe matter in a strong 
_B-field. As a standard reference for our comparison we 
use the modified Urea process. The presence of a magnetic 
field has two different effects on neutrino emissivities: (i) 
in pure neutron matter it allows for spin-flip transitions 
where the finite difference of the momenta of neutrons at 
two different Fermi surfaces enables one to satisfy energy- 
momentum conservation, (ii) the charged particles occupy 
Landau levels leading to a smearing of the transverse mo- 
menta over an amount \/e_Bp]. As a consequence there are 
two typical scales of the magnetic field, where effects on 
the emissivity are expected; first for |/is|-B ~ T which is 
relevant in neutrino-pair bremsstrahlung from neutrons, 
and second eB ~ p 2 F relevant for the Urea process. 

As well know, the emissivity of any particular reaction 
can be related to the imaginary part of the polarization 



function of the medium (Voskresensky & Senat orov 1 1986 , 
Raffelt & Seckel |l995| , Sedrakian & Dieperink |1999| ). We 
compute the polarization function of neutrons and protons 
in strong magnetic fields employing the finite temperature 
Matsubara Green's functions technique. For the case of the 
bremsstrahlung the time- like properties of the polarization 
function are relevant. The space- like properties of the po- 
larization function, relevant for the ncutrino-nucleon scat- 



tering, have been studied by Arras & Lai ( 1999 ) in an 
equivalent response function formalism. 

The plan of this paper is as follows. The 
bremsstrahlung emissivity is related to the polariza- 
tion function of the medium in Sect. 2. The neutrino 
emissivity via neutrino-pair bremsstrahlung from neu- 
trons is discussed in Sect. 3 and that from protons in Sect. 
4. The Urea process in strong magnetic fields is briefly 
discussed Sect. 5. The EoS of npe matter in a magnetic 
field is discussed in Sect. 6. Our numerical results are 
presented in Sect. 7. Sect. 8 contains our conclusions. 



2. The bremsstrahlung emissivity 

The neutrino emissivity of an infinite medium of inter- 
acting hadrons can be expressed in terms of the imagi- 
nary part of the finite temperature polarization function 
n(q, to), where q and u> denote the momentum and energy 
transfer to the leptons. In the single-loop approximation 
the finite temperature polarization function in a magnetic 



field is a 2 x 2 matrix in spin space (Mattuck |1992| ) 
-2 f d 3 p 

i 

i(p+q) 



—2 f dp 

IL 3 , a ,(q,u) = — J ^—3 G s ,(p + q,ip + iq) 

G s (p,ip), (1) 



where (3 is the inverse temperature. The single-particle 
Green function can be expressed as G s {p 1 ip) = [ip — 
(E PtS — Ep)}^ 1 , where Ep is the Fermi energy and ip is the 
Matsubara frequency; s',s = ±1 specify the initial and fi- 
nal nucleon spins. (We assume the magnetic field along the 
positive z— axis.) Carrying out the frequency summation 
and taking the imaginary part one finds the well-known 
result (Mattuck |1992[ ) 



ImlL,^ (q, w) = 2n J 



d 3 p 



(27T 



3 (f( E p,s) - f{Ep+q,s')) 



3{Ep, s — E, 



p+q,s' 



(2) 



where f(E PtS ) = [ e 0(Ep,s-Er) + is the Fermi- 

distribution function. The emissivity is then given by 

^ = 3 £ |M S , S , | 2 / §^§^d. A ImII s>s , (g, U ) 

ss' 

g B {uj)u)5(u-uj l ,-uj 17 )5 3 (q-p l ,-p 17 ), (3) 

where gs (w) = [e ~ 1] _1 is the Bose distribution func- 
tion, the factor 3 comes from the sum over the neutrino 
flavours, and the weak interaction matrix (neglecting lep- 
ton momenta) is | Af SjS / 1 2 = (G F /2) 2 (c v S SjS > + c A (25 s .- s ' + 
$s.s')) with cy and ca the vector and axial vector coupling 
constants and Gf the Fermi weak coupling constant (see 
Appendix B). 



3. Neutrino-pair bremsstrahlung from neutrons 
(n — ► n + v + 17) 

In the absence of a magnetic field the imaginary part of 
the polarization function vanishes for time-like processes 
in the quasi-particle approximation, because energy and 
momentum cannot be conserved simultaneously. In a pure 
neutron system a magnetic field B will give rise to a dif- 
ference between the Fermi momenta of the neutrons with 
spin parallel and spin anti-parallel to the i?-field (see Ap- 
pendix A) 



PFn 



(E 



sg„eB) 



1/2 



We use the natural units, h 



where u n — g n e/(2m) is the neutron magnetic moment 
and g n = —1.91 is the neutron g- factor and we assume 
fi n B <C m. For fi / energy-momentum conservation 
can be satisfied, and as a result one may expect that a 
field with strength |/x„|i? ~ T leads to a finite spin- flip 
polarization function whenever u> ~ 
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3.1. Polarization junction 

In the time-like region it is preferable to use the relativis- 
tic kinematics. The non-relativistic kinematics does not 
produce the correct zero-field limit H(q 7 uj,B = 0) = 0, 
rather a spurious finite contribution H(q,ui, B — 0) oc 
exp(— m/T). To evaluate the angular integral in eq. (||) 
note that the energy conserving ^-function is non-zero if 
| cos{6 pq )\ < 1, where 9 pq is the angle between the momen- 
tum vector p of the neutron and the momentum transfer 
vector q; this condition yields the minimal and maximal 
values of the three-momentum of neutrons (p m in and p max ) 
for which the imaginary part of the polarization function 
is finite. 

To obtain a real solution for p m i n and p max using the 
relativistic energy/momentum relation, E 2 s = m 2 + p 2 — 
sg n eB, the conditions a = [(4m^, n B + u> — q 2 ) 2 — 4(m 2 + 
2mfi n B)(uj 2 -q 2 )] > and s' = -s = 1 should be fulfilled. 
The result is 



Pmin(g,w) 



and 



-4m/i„B — lu 2 + q 2 )q — LOy/a\ 
2{oj 2 -q 2 ) : 



-4mfi n B — u) 2 +q 2 \q J r LJ^/a 
2(uj 2 -q 2 ) 
00 



b > 
b < 



(4) 



(5) 



with b = 2Ep jS uj 



AmiJL n B + {J 2 — q 2 (note that \i n < 
0). Replacing the integral over the absolute value of the 
momentum p by an energy integral dE/dp = p/E w p/m 
yields 

ImII s>s /(qf,w) = — / -z—(f(Ep,s 
q .L . 2n 




with E„ 



^mm/max ( m T ^ m i n / ma x 

For large j3 the rhs of Eq. (0) 



L n 



(G) 



- sg n Be)^ 2 - E Fn . 
is non-negligible only if 



-Emm < and E max > 0. The latter condition is of minor 
importance, since it is satisfied for almost all q and u>. 
As is illustrated in Fig. |l| for small T the region in the q, ui- 
plane where qlmll(<7, uS) is finite, is essentially bounded by 
three straight lines (which become exact limits for T — > 0); 
these boundaries are essentially determined by the fact 
that q < to and the condition p m m < p Fn ■ The latter 
can also be expressed (neglecting all terms in p m - ln except 
the leading order terms in m) as u> > h~ ,uj < h + with 
ft* « (2mW\B±p- Fn q)/(m 2 + fej 2 ) 1 / 2 . 

3.2. Emissivity 

We have only to consider the case s' = — s = 1, in which 
case M 2 S , = G 2 F c 2 A /2. By integrating over the neutrino 
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Fig. 1. Contour plot of qhaIl StS i (q, uS) in units of (m 2 /%) 
for s' = -s = 1 at B = 10 16 G, T = 10 9 K and density 
n = 0.155 fm -3 with x = (3q and y = 0w. 



momenta in Eq. (g) and introducing the dimensionless 
parameters y = (3lu and x = f3q, the emissivity can be 
expressed as 



n2 2 2 poo 4 



-0E mln -y + i 



1 -m ' 







1 



/3B„ 



1 



(7) 



To obtain some insight in the dependence of the emissiv- 
ity on B we distinguish three different regions of B (for 
simplicity we take \q\ = to, so that the integral over x can 
be replaced by 2y/5) 

— For weak magnetic fields (\nn\B <C T) the region in 
which E m i n < is proportional to /i n B/T and peaks 
around y w \i n B jT '. For y <C 1 one obtains 



G F c\m 2 7 6 
30(2^) 5 ^ ma * 



Vmin ) ) 



(8) 



where y ma x/ mm = [2fj, n B/mT) {^J m 2 + (p Fn ) 2 ± 
p Fn ). Hence for fixed T, one finds 

e,T7 oc B e . (9) 

— For strong magnetic fields (|/z n |i? ~ T) the Bose func- 
tion in Eq. must be kept and the emissivity be- 
comes 



J F A 



5(2tt) 5 



-T 



dy 



y 



ev - 1 



(10) 
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The integral peaks around |/i n |-B ~ 3T with a width 
6T(p Fn /m). As discussed in section 7 for B values in 
the range \fi n \B ~ T the emissivity becomes compara- 
ble to conventional (modified Urea) process. 
— From Eq. ^ we sec that for superstrong magnetic 
fields (\nn\B S> T) the emissivity falls off exponentially 
with the magnetic field for fixed T. 

4. Neutrino-pair bremsstrahlung from protons 
(p — ► P + v + V) 

In a magnetized matter in addition to the spin-magnetic 
field interaction the charged particles (protons, electrons) 
are grouped into Landau levels. The proton's Fermi mo- 
mentum in a magnetic field for given Landau level number 
N and spin s is given by (see Appendix A) 



(E 2 Fp -m 2 -(2N+l-sg p )eB) 1 / 2 , 



s,N 

Pf p 

where g p is the proton g-factor. 

The population of Landau levels leads to a smearing of 
the transverse momentum, Apj_ ~ VeB, and as a conse- 
quence the condition for energy/momentum conservation 
is softened. In the special case of a superstrong magnetic 
field, eB > (p Fp ) 2 , only the lowest (N = 0) Landau level 
is occupied. 

4-1- Emissivity 

In the present case it is convenient to define a reduced 
proton polarization function for a specific Landau level 
N = p\/{2eB) 

lmU N ' >NiS ^ s (q z ,UJ z ) = J dp z [f(E pZ:N:S ) - 

f(E pz ',N',s')] ${E pz ^N,s — Ep Z ^N'.s' + LU Z + Un',n), (11) 

where lo z = (u> 2 - i^ 2 _) 1 ^ 2 with lu± = (q 2 + q 2 ) x l 2 , Un'.n 
is the energy difference between the Landau levels N, N' 
and the proton energies E pz ^N,s are defined in Appendix 
A. In evaluating Eq. ( fil| ) the integral over p z is replaced 
by (m/ppp) J dE pz N,s, and after integration over E pz ^ s 
one obtains 



lmil N ^ N ^i. s (q z ,uj z ) 



— lf(E a ) + f(E b ) 
Pf p 

-f(E a + uj z + U N , N ) - f(E b + uj z + U N , N )] 
with 

E a ,b = (m 2 + pl b + (2N + 1- sg P )eB) 1 l 2 , 



(12) 



(13) 



where p a and p b follow from the condition of conservation 
of z-momentum and energy. To obtain real solutions for p a 
and pb the conditions d = (q 2 — u 2 + AmpLpB) 2 — Am 2 (uj z — 
q 2 ) > (neglecting some small terms) and s = — s' = 1 
should be satisfied. The result for p a and pb is 



Pafi(q Zl LJ z ) 



(AmppB + q 2 - uj 2 )q z ± u z Vd\ 



(14) 



We note that the polarization function essentially vanishes 
unless p a and pb are close to the Fermi momentum p Fp — 



(PFp) 2 ~ 2NeB. In the case of a superstrong magnetic 
field only the N' = N = states are occupied and the 
contribution to lmH^i^^ s ^ s (q z ,uj z ) comes only from the 
lines in the w z ,g z -plane defined by (neglecting all terms in 
p a except the leading order terms in m) 



2m/j, p B ± q z pp° p 



2(1 - g p )mn B B + (p F ° p ) 2 



(15) 



For weaker S-fields a larger space in the q z ,w z -plane con- 
tributes eventually leading to a situation similar to that 
for neutrino-pair bremsstrahlung from neutrons. 
The emissivity can be written as 



^ 4(2tt) 

N',N,s',s y ' 



^-^ J divdq z duj z g B {uj) 



lmR N > <N<s > jS (q z , u z ) \M N >. NiS i tS \ 2 , 

where 

\M N , tNiS ,, s \ 2 =J2\ M R',R, 



(16) 



,N'N,s',s\ 



R'R 



eB 
2^ 



2c A S s >- s + c 2 A 6 s > tS ]lff, N (v) 



r 2 

with v = (lu 2 — oj 2 )/(2eB), R and R' the guiding cen- 
ter quantum numbers and In,n' the associated Laguerre 
polynomials (see Appendix A). As in the case of neu- 
tron bremsstrahlung only one spin configuration con- 
tributes (s = — s' = 1). For strong magnetic fields, 
B ^> T 2 1 (4to/xb), the dominant contribution to the emis- 
sivity comes from N' = N in which case the weak inter- 
action matrix can be simplified in case of Nv <C 1 



E 

R'R 



\-^R',R,N,N,- 



\ 2 -9}lr 2 —r 2 (v) 



G 2 F 2 eB 



2 Ca 2^- 



(17) 



Using the fact that Un'n vanishes for N' = N, in this 
case the emissivity becomes 



_ ^ 2Amc 2 A G 2 Fi i B BT 7 f 1 y 4 

' vv ~ 2^ 4(2tt)5 J y e y-i 6 

N=0 v ; J 

y ryz 

dy z / dx z lmH N:N _ t+ (x z ,y z ), 
o Jo 



(18) 



where y = u/T, y z = uj z /T, x z = q z /T and N max = 
(p~Fp) I '(^ m l J 'BB) . In weak magnetic fields when eq. ( |i"7j ) is 
not a good approximation, the summation over N must 
be restricted to (4m/xsB)/ '(36T 2 ). In general the integral 
over x z must be carried out numerically. In case of fi p B « 
T the main contribution comes from y z — y, and as a 
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result the emissivity of neutrino-pair bremsstrahlung from 
protons is very similar to neutrino-pair bremsstrahlung 
from neutrons. On the other hand for weaker fields the 
emissivity of neutrino-pair bremsstrahlung from protons 
is larger due to the smearing of the transverse momenta 
of the protons. 

5. Direct Urea process (n — ► p + e + V e ; 
p + e — ► n + v e ) 

It is well known that the direct Urea process can occur 
only if the Fermi momenta satisfy the inequality 



PFp + PFe > PFr. 



(19) 



(except for a small thermal smearing), i.e. the proton con- 
centration x p — n v / (n v + n n ) needs to be larger than 1/9 
(L at ti mer e t al. 1991 ). Recentl y it w as shown (Leinson & 
Perez |l99S| ; Baiko & Yakovlev |l999| ) that the presence of 
a magnetic field has an effect on the Urea process; in par- 
ticular it leads to a non-vanishing emission rate in case 
the triangular condition ([H]) is not satisfied (the so-called 
classically forbidden region). This is caused by a smearing 
of the momenta in the presence of a magnetic field which 
leads to a softening of the sharp border between the al- 
lowed and forbidden regimes. Since the emissivities due 
to th e Urea process in superstrong 5-fields (Leinson & 
Perez 1998) and in arbitrary B- fields (Baiko & Yakovlev 



199G) have been derived previously, we provide below a 



brief discussion for the purpose of completeness. 
5.1. Emissivity 

The emissivity is calculated using the quasi-particle ap- 
proximation for the polarization function where the elec- 
trons and protons occupy the Landau levels N e , N p , re- 
spectively, 

e vV = 2 ^ / <^Pnd,Pez<^Pvdp V z 

N a N p s n s p ' n > 

${Pnz ~ Ppz ~ Pez ~ Pvz) ^ | M Rei R pi N e ,N p ,s n ,s p | 2 
R S R P 

w v f(E n )(l-f(E?))(l-f(E e )) 

5{E n + V np -EP-E e + Lj u ), (20) 

where V np — E Fp + E Fe — E Fn and the neutron, proton, 
and electron energies are given in Appendix A and the 
matrix element M. in Appendix B. The factor 2 comes 
from the contribution of the inverse reaction. 

5.1.1. Superstrong magnetic fields 

First we consider superstrong magnetic fields where the 
electrons and protons populate the ground state Lan- 



that \p pz \ w p F ° p and \p ez \ w p° Fe . Neglecting the neu- 
trino momentum the matrix element in Eq. ( p0| ) simpli- 
fies to T,r c r p \M\ 2 = Cs n cxp(-plJ2eB) with C s „ = 
G 2 F [(l-2cA + 4)<5 Sn , 1 + 4 C ^ s „,_ 1 ]/2,andp2 ± = (pJ-J*- 



(Pp 



Evaluating the energy integrals leads to 



= > a 



4(120 



(eO/ F «„)exp[-(^) 2 /2eB] 



(2tt)5 



6w 2 )Be m n m p ^ 6 



D +o 

PFp 



\p 



Fp ' 



PFe 



«P[-((P&) 2 " (p + F ° p +P Fe ) 2 )/2eB}). (21) 

The Q functions correspond to conservation of momentum 
in the z-direction; it is worth noting that the triangular 
condition expressing the momentum conservation in the 
presence of a superstrong B- field, \p Fp + p° Fe \ < \p s Fn \, is 
the opposite of the one found for the B=0 case. 
The expression (^) for the emissivity has been ob- 



tained previously by Leinson & Perez (199S) and Baiko 



k Yakovlev ( p9S| ) . The factor C+ = GplT - c A ) 2 /2 i n 



Eq. (^Tj) agrees with the result of Baiko & Yakovlev (1999) 
but differs from the one given by Leinson & Perez (1998) 



5.1.2. Weak magnetic fields 

In the case of weak magnetic fields, i.e. eB < p Fp , the 
summation over the Landau levels in Eq. ( pp| ) is replaced 
by an integral^. Also the Inn functions in the matrix el- 
ement given by Eq. ( |B.9| ) of Appendix B are replaced by 
their small- -B asymptotes 



2eB T kN p (v) 



N C N P 



-^Q e J d P 2 e±dp 2 p ±A{p p ±,p e ±,B)Ai 2 (z) 
with Ai(z) being the Airy function with 

_ [Pl± ~ (Pp± +Pe±) 2 ](Pp± +Pe±) 1/3 



(22) 



v = plj_/2eB, and 

A(Pp-LiPe±, B ) = 



2 e S) 2 /3(p p±+J3e± )4/3 

(2eB) 2 / 3 



(Pp± + Pe±) 2/3 (Pp±Pe±) 1/3 ' 

Neglecting the neutrino momentum in the z-direction 
in the delta function in Eq. (|2p| ) in comparison with the 
momenta in the z-direction of the other particles and per- 
forming the integral over d co s(0„), the emissivity can be 
written as (Baiko & Yakovlev |1999 ) 

8(1 + 34)G F (120 + 6ir 2 )m n m pPFp p Fe 
tvv ~ (2TrfeB 

J dcos(9 p ) dcos(9 e ) A(p p j_,p e j_, B) Ai 2 (z) 

5{pFn - \pFpCOs{9 p ) + p Fe COs(9 e )\). 



(23) 



dau levels (Af e 



N„ 



0) with spin parallel to B, so 



2 In this Sect, we use the definition of the Fermi momenta 
as in the field free case. 
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In the classically forbidden domain for the direct Urea pro- 
cess for finite B the reaction becomes possible due to the 
tunnelling mechanism. In the quasi-classical approach the 



emissivity can be expressed as (Baiko & Yakovlev 1999) 



forbidden. 

e — 

vv 



with 

R(x,y) 



(B?0)=R(x,y) e%° wed (B = 0) 



x + 12y x 3 / 2 



exp(- 



,3/2 



-). 



(24) 



(25) 



(26) 



where y = N F ^ = (p Fp / (2eB)) 2 / 3 and x is a measure of 
the violation of the triangular condition, 

_ P% ~ (PFp+PFe) 2 

X ~ v 2 N~ 2/3 ' 

PFp^Fp 

One sees that with increasing x, the emissivity decreases 
exponentially, and therefore the effect is important only if 
PFn does not exceed (jpFp +PFe) significantly. 
On the other hand in the classically allowed region, where 
the inequality is satisfied, the -B-field will give rise to small 
quantum oscillations of the emissivity. 



6. Equation of state (EoS) 

In order to estimate the effect of a magnetic field on the 
various cooling processes in a neutron star and to compare 
with the conventional result we employ a simple model 
EoS. It is assumed that the neutron star matter consists 
of neutrons, protons and electrons only (npe-matter) with 
two conditions imposed: charge neutrality, n p = n e , and 
/3-equilibrium, </>„ = cf> p + (fijj^ The EoS consists of the en- 
ergy density as a function of the densities of the particles. 
The non-relativistic energy density can be decomposed as 



(Lattimer & Swesty 1991; Balberg & Gal 1997) 



E — Ef-in + E mag + E mass + E pot + Ei e p. (27) 
Here the kinetic energy density is the sum of the neutron 



pn _ 3 ^ (Pirn) „a„ 
km ~ 5 ^ 2m„ " 

and proton contribution 

P _ly (Pfp°) 2 Sb 
km ~ 3 f" 2m p p ' 

where the proton density for spin s p is 

s _ \ - eB_ 

lpP ~ ^ 2ir 2 

N 



//?.'■ = > — \Hp h ff) 2 -2NeB, 



(28) 



(29) 



(30) 



because for finite B the integrals over the transverse mo- 
mentum of the proton with respect to the magnetic field 

3 <j>i is the chemical potential of a particle; it can be obtained 
from the relation <f>i = dE/drii. 



can be replaced by a summation over the Landau levels. 
The interaction of the magnetic field with the spin of the 
neutron and the proton is 



E, 



mag 



n B Sn<" - n P B Y Spn p " 



(31) 



The E mass term contains the masses of the two nucleons 
E m ass = m n 2J n s ™ +m p ^2 n P p ■ ( 32 ) 

The Ep t term is the potential energy density, which is 
parametrized as (Lattimer & Swesty 1991) 



Epot = an 



bn 



l+d 



(33) 



where n — n p + n n is the total nucleon density. The last 
term on the rhs of Eq. (p3|), which corresponds to the 
symmetry energy, influences to a large extent the proton 
fraction n p /(n p + n n ). Finally the energy density of the 
electrons is 



Eiep — Ep e n e , 
where 

E9N c eB 



(p° Fe ) 2 - 2N e eB 



(34) 
(35) 



with gjv e = 1 for N e = 0, g^ e = 2 for N e > and N e is 
limited from above by N™ ax = {p° Fe ) 2 / (2eB). 
We obtain the EoS of the npe-matter in magnetic fields in 
the standard manner, by assuming a given nucleon density 
n and T — and solving the equations of charge neutral- 
ity and /3-equilibrium. In two limiting cases the solutions 
are straightforward. 

i) B = 0. The summations over the proton and electron 
concentration in Eqs. (^) and ( |35| ) can be replaced by 
integrations and the densities n e , n p and n n can be cal- 
culated in the standard way using m = p Fi /(3ir 2 ). Using 
the /3-equilibrium condition, 4> e = <f) n — (f> p , the chemical 
potential of the electron is given by 



(3^ 2 n„) 2 / 3 (3ir 2 n p ) 2 / 3 



2m r , 



2m»- 



+ 4c(n p - n n ) + A (36) 



with A the mass difference between neutron and proton, 
ii) Superstrong magnetic fields, (p F p) 2 < (2eB), in which 
case the protons and electrons are in the ground state 
Landau level. Then the electron chemical potential is 



2 n M2/3 



1 ^ (6tH<") 
, 9p ~ 1 



1 /2TT 2 n r ,\ 2 



2m r 



\ eB ) 



2m r 



eB + Ac(n p — n Tl 



(37) 



In this case the proton fraction depends on the magnetic 
field. In Eq. @ we take d = 2, a = -285.1 MeV fm 3 , 

-107.1 MeV fm 3 (Lattimer & 



b = 968 MeV fm 6 , c 
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from protons 
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Fig. 2. Proton fraction as a function of density for various 
B-ficlds 



Fig. 3. Emissivities of various processes at saturation den- 
sity n — no and for B = 10 16 G. The contribution of the 
direct Urea process is negligible. 



Swesty |l99l|) . 

In Fig. |2| the resulting proton fraction in the absence 
and in the presence of magnetic fields is plotted. The 
kinks in the proton fraction at certain densities corre- 
spond to the occupation of a next Landau level (see also 
Suh & Mathews 1999). For superstrong magnetic fields 
these kinks strongly influence the proton fraction, but for 
weaker magnetic fields (B < 10 17 G and n > no where 
no = 0.155 fm~ 3 is the saturation density) they do not 
affect the proton fraction. Our results are in agreement 
with the previous results for the proton fraction derived 



from different equations of state (Lai & Shapiro 1991; Suh 
& Mathews ^999] Broderick, Prakash and Lattimer |2000| ) . 



7. Results 

The emissivities for the various energy loss processes are 
compared in Figs. [| || and || for n = no for three differ- 
ent magnetic field strengths (10 16 ,10 17 and 2 x 10 18 G, 
respectively). In the first two values of the B-field the 
EoS based on B — is used (because in these cases the 
influence of the magnetic field on the EoS itself is small) ; 
for the third value of the -B-field we use the EoS for su- 
perstrong magnetic fields. 

To enable a comparison with the neutrino processes which 
are routinely included in the cooling simulations of neu- 
tron stars, we show also the emissivity due to the modi- 
fied Urea process in the zero field limit. The relevant ex- 
pression for the emissivity of the modified Urea process, 
without corrections for the magnetic field, and valid for 
m* — m is (Friman & Maxwell 1979) 



so i-io — 

o 



in 

I i-io 24 h- 



— i r 

Direct Urea 
. Modified Urea 

Bremsstrahlung 
from neutrons 
Bremsstrahlung 
from protons / / / 




= 1.8 x 10 21 (n/n Q ) 2/3 T { 



g ergs cm 3 s 1 



(38) 



T(K) 



Fig. 4. Same as in Fig. 3 but for B = 10 17 G. 



with Tg the temperature in units of 10 9 K. 

The temperature region where the one-body neutrino- 
pair bremsstrahlung is important increases with increasing 
mag netic field (Figs. | and |). The pair bremsstrahlung 
from neutrons is efficient whenever |/Lt n |-B ~ T, since then 
the energy involved in the spin-flip is of the same order of 
magnitude as the thermal smearing of the Fermi surface. 
The temperature at which neutrino-pair bremsstrahlung 
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Fig. 6. Emissivity at n = 3n and B = 10 lb G 



from neutrons becomes comparable to the competing pro- 
cesses roughly coincides with this condition. For lower 
temperatures the emissivity drops exponentially, because 
the energy transfer becomes larger than the thermal 
smearing. Neutrino-pair bremsstrahlung from protons is 
important when [i v B ~ T . The emissivity due to the pro- 
tons increases faster than the emissivity due to the neu- 
trons with the temperature, as the smearing of the pro- 
ton transverse momenta provides an additional relaxation 
on the kinematical constrains. However for temperatures 
smaller than [i v B ~ T the emissivity drops just as for 
neutrons exponentially. 

As seen from Figs. || and |] the emissivity of the mod- 
ified Urea process is larger than that of neutrino-pair 
bremsstrahlung from neutrons and protons at high tem- 
peratures mainly due to the different temperature de- 
pendencies of these processes (cx T 7 for the one-body 
bremsstrahlung as compared to oc T 8 for the modified 
Urea) . In the case of a superstrong magnetic field the large 
uncertainty in the transverse momenta of the protons and 
electrons allows the direct Urea process to occur (see Fig. 
|^) and its emissivity dominates the emissivity of any other 
processes. In Fig. |^ the emissivities are shown for n = 3no 
and B = 10 16 G. At this density the triangle condition 
Pf p + PFe > PFn is satisfied, so that the direct Urea pro- 
cess is allowed and dominates the cooling. The emissivities 
of the other processes are slightly larger than those shown 
in Fig. U due to the fact that the density is larger. 

8. Conclusion 

We have studied the neutrino emissivity of strongly mag- 
netized neutron stars due to the one-body processes driven 



by the charged and neutral current couplings between 
the neutrinos and baryons. We have shown that, in ad- 
dition to the well-known charged current process (the di- 
rect Urea reaction), there is a new channel of energy loss 
- the one-body neutrino pair-bremsstrahlung in a mag- 
netic field. The process does not have an analogue in the 
zero field limit and competes with the modified (two-body) 
bremsstrahlung process as the dominant neutral current 
reaction for fields on the order 10 16 — 10 17 G and temper- 
atures a few times 10 9 K. For superstrong magnetic fields 
in excess of 10 18 G the direct Urea process takes over. 
Our numerical evaluation of the emissivities of sev- 
eral competing reactions, which is based on a simple 
parametrization of the EOS of npe-matter in a strong 
magnetic field, shows that under certain conditions the 
emissivities of the one-nucleon processes, such as the di- 
rect Urea and the one-body bremsstrahlung, are of the 
same order of magnitude or dominate the standard pro- 
cesses commonly included in the cooling simulations in the 
zero-field limit. 



Appendix A: Electrons, protons and neutrons in 
a magnetic field 

The wave equation for a fermion with charge q and mass 
m can be written as (Itzykson & Zuber 198C ) : 



(i@ + q4-(m- ^-iib^ v F^))^[x) = 0, (A.l) 

with fis = e/2m, /6 — 7 M <9 M , A = 7 M A' i and Ag the 
anomalous gyromagnetic factor. Here 7 M is a Dirac matrix, 
A M the vector potential, Fa V = d^A v — d u A^ and a^ v = 
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r[7 M ,7"]- Or equivalcntly 



A. 2. Proton 



{{id + qA) 2 + 1cj^F^- 

In particular, if one chooses the Landau gauge A^ = 
(0, — ^p, ^p,0) with the magnetic field in the ^-direction, 
one finds 

((id + qA) 2 + sqB - (m - Agsp B B) 2 )^{x) = 0. (A.3) 



A.l. Electron 

For an electron (q — — e and Ag = 0) the energy eigenval- 
ues are given by 



E N ,„ = (m 2 +pl + 2NeB) 1 ' 2 , 



(A.4) 



with N denoting the Landau level number. The eigenfunc- 
tions are factorized (Sokolov & Ternov 1968; Arras & Lai 



1999| ) 



^ = e i **- iE »>- t Ue{ P ,, 



iVblN,R.(t) 

V iV d I N , R (t) ) 



(A.5) 



and 



m 



ww = vi e ~ t/2 ' M/2£ ^(*)' 

where L% (t) 

_ j L%(t) if M > 

-{(-l)mt^\L^ m (t) if M<0 

with t = p 2 eB/2, V a = C+D+, V b = aC-D+, V c 
o-C+D-, V d = C-D_, 



c = _L( 1±(7 p* 

± V2\ {p 2 z +Am^ B BNy/ 2 



1/2 



D+ = -4= I i ± ^ 



and cr is the longitudinal spin projection along p + eA; 
L|f (t) is the generalized Laguerre function. The Landau 
level number N and the guiding center quantum number 
R are positive integers. (Note that in the lowest energy 
state (N = 0) the spin s can only have the value -1). 
In matter in the presence of a magnetic field the Fermi 
energy is given by E F = (to 2 + {Ppe) 2 + 2NeB) 1 ^ 2 so that 
a magnetic field will give rise to a difference in the Fermi 
momenta pp e . The degeneracy is 1 in case of N = and 
2 in case of N > 0. 



For a proton (q = e and Ag = g p — 1 with g p = 2.79) the 
energy eigenvalues (neglecting the i3 2 -term) are 



( A - 2 ) E NjS = (to 2 + p\ + (2N + 1 - sg p )eB) 1 ' 2 . 



(A.6) 



The non-relativistic eigenfunction of a proton in a mag- 
netic field is 



ipzZ-iE N i3 t 



U p (p, 4>) = 



U P (p, (f>), 



J(R-N)<j> 



{Ss,llR,N(t) \ 
5 s -llR,N(t) 





(A.7) 



(Note that the role of R and TV are interchanged compared 
to the electron case.) In matter the Fermi energy is given 
by E F = (to 2 + (p^) 2 + (2N + 1- sg p )eBf/ 2 . 



A.3. Neutron 

For a neutron (q = and Ag = g n with g n = —1.91) in 
not extremely large magnetic fields the B 2 -term can be 
neglected; then the energy eigenvalues are 



E s = (to 2 + p 2 - 2sg n mpBB) 



1/2 



(A.- 



and the eigenfunction are simple plane waves. In neutron 
matter in thermal equilibrium the Fermi energy is given by 
Ef = (m 2 + (ppn) 2 ~ 2sg n mpBB) 1 / 2 . Hence the neutrons 
with spin up (down) occupy two Fermi spheres with Fermi 
momenta related by (p% n ) 2 = (p Fn ) 2 + 4m/x n S. 

Appendix B: The weak interaction matrix 
elements 

B.l. Neutrino-neutron interaction (neutral current) 

The weak interaction matrix for the neutron 
bremsstrahlung is 

G F 



AL 



2V2 



dr V'„i7^(cy-CA75)^n2V' 1 ,7 AI ( 1 -75)'0F-( B - 1 ) 



Here cy and ca are vector and axial-vector coupling con- 
stants. The interaction matrix consists of a hadronic part 
and a leptonic part 



\M sy \ 2 = ^X^L^ 

with the leptonic part given by 

U»> = JL [py* _ g ^(p . p>) + p'pin + ie^paj/p] 
pp' H 

and the hadronic part by 

I c^Ss^ 

c 2 A 5 s - s , 

— icr A s8 St - s i 

\ -s'cvc A 5 s , s i 



(B.2) 



Xuu — 











-s'c V C A 5s,s' \ 




c A 5 Si s' J 
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Contracting the hadronic and leptonic tensors and ne- 
glecting terms, which after integration over the neutrino 
momenta vanish one finds 



Vi (i — a,b,c and d) can be approximated by D « \j\[2. 
This leads to the following expression for the leptonic part 



with Cy 



r 2 



c 2 A )8 s 



2c 2 A 6 3 , 



(B.3) V v 



1 and ca = 1.26. 



( J+{v) J-(v)\ 

J+(v) iJ_(v) 

-U-(v) J + (v) 

\J-(v) J+(v)J 



(B.6 



-B.,2. Neutrino-proton interaction (neutral current) 

In case of protons the integration over space coordinates 
requires special attention, because the wave functions of 



with J±(v) — /jv e ,iV ( w ) ± In c -i.n p ( v )- As a res ult the 
following relation is obtained for the squared norm of the 
interaction matrix 



protons are not simple plane waves (Eq. A. 7). Defining 
u± = (Pux + Pvx)e x + ip vy +Pr7y)e y , x ± = xe x + ye y , 
v = uj 2 _/2eB and carrying out the summations over R' 
and R, leads to 



\M Rp , Re ,N p . 



N, 



7? e H, v 



CPeB 
~2~1^ 



(^ P , S „(cy + c\) 



I N e -l,NS V ')) 



+28, 



- ■ s n 9a,Sp(Ijy e -i,N, 



N C ,N, 



.(«)) 



R'R 



R'R,N',N,s's\ 



2tt 



R',R 

pdpIw,N>{t)I R ,N(t) 



eB 

\M s ', s \ 2 —I 2 N ,, N (v) 
Zir 



- 2 S Sp ,- Sn gl(6 Sp>1 I 2 Ne M{v) 



^.- ltJV »))(B.9) 



with cy = 1 and ca = 1.26. 



(B.4) 



with cy = —0.08 and ca 



-1.26. 



B.3. Direct Urea process (charged current) 
The interaction matrix is 



M 



R e ,R p ,N e ,N p ,s n ,s p 



Gp 
V2 



dripp^^Cv - CA"f5)lpn 

^(1-75)^-. (B.5) 



The eigenfunctions of the electron, proton and neutron 
are given m Appendix A. The function Ir N p (t) with 
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